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Abstract 

The complexity of the following numerical problem is studied in 
the quantum model of computation: Consider a general elliptic partial 
differential equation of order 2m in a smooth, bounded domain Q C M d 
with smooth coeficients and homogeneous boundary conditions. We 
seek to approximate the solution on a smooth submanifold M Ç Q of 
dimension < d\ < d. With the right hand side belonging to C r (Q), 
and the error being measured in the L 00 (M) norm, we prové that the 
n-th minimal quantum error is (up to logarithmic factors) of order 

- min((r+2m)/di, r/d+l) 

For comparison, in the classical deterministic setting the n-th minimal 
error is known to be of order n~ r ^ d , for all d\, while in the classical 
randomized setting it is (up to logarithmic factors) 

- min((r+2m)/di,r/d+l/2) 



1 Introduction 



The complexity of solving elliptic problems in the classical deterministic 
setting was studied in [30, 31, 7, 5, 6]. In [18] such problems were considered 
in the classical randomized setting. The quantum complexity of ordinary 
differential equation was investigated in [19], while in [21] certain parabòlic 
problems were studied in this setting. The complexity of elliptic problems in 
the quantum model of computation has not been analyzed before. This is the 
tòpic of the present paper. We consider a general elliptic partial differential 
equation given on a smooth domain in W d , with smooth coefflcients and 
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homogeneous boundary conditions. We seek to find an approximation to the 
solution on a given, cíi-dimensional smooth submanifold, where < d\ < d. 
Thus, we consider the whole range of problems from local solution (find 
the solution in a single point, d\ = 0) up to global solution (find the full 
solution, in the whole domain, d\ = d). Our analysis is carried out in the 
quantum setting of information-based complexity theory, as developed in 
[11]. For a study of other bàsic numerical problems in this framework we 
refer to [23, 12, 14, 15, 19, 21, 27, 32], see also the surveys [13, 16]. For 
general background on quantum computation we refer to the surveys [2], 
[8], [26], and the monographs [25], [9], [24]. For the classical settings of 
information-based complexity theory we refer to [28, 22, 10]. 

This paper can be considered as a continuation of [17, 18]. The approxi- 
mation of weakly singular integral operators plays a key role again. In some 
situations, techniques from [17, 18] can also be applied to the quantum set- 
ting, while in others entirely different approaches are needed. In particular, 
a number of new tools for the general quantum setting of information-based 
complexity has to be developed. 

The paper is organized as follows. In section 2 we describe the quantum 
setting, general results about quantum n-th minimal errors are derived in 
section 3. In section 4 we study weighted mean computation and integra- 
tion. These are preparations for section 5, in which we are concerned with 
quantum approximation of weakly singular operators. Section 6 contains 
the statement and the proof of the main result about the query complexity 
of elliptic PDE. Finally, in section 7 we recali the respective results of the 
classical deterministic and randomized settings and compare them with the 
quantum setting. 

2 Notation 

A numerical problem is given by a tuple V = (F, G, S, K, A) , where F is a 
non-empty set, G a normed space over K, where K stands for the set of real 
or complex numbers, S a mapping from F to G, K a non-empty set and A 
a non-empty set of mappings from F to K. We seek to approximate S(f) 
for / e F by means of quantum computations. 

Usually F is a set in a function space, S is the solution operator, which 
maps the input / £ F to the exact solution S(f), and we want to approx- 
imate S(f). A is usually a set of linear functionals, supplying information 
A(/) about / through which the algorithm can access the input /. K is 
mostly K or C, G is a space containing both the solutions and the approxi- 
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mations, and the error is measured in the norm of G. 

In the sequel it will be convenient to consider / G F also as a function 
on A with vàlues in K by setting /(A) := A(/). Let ^"(A, K) denote the set 
of all functions from A to K. 

Let Hi be the two-dimensional complex Hilbert space C 2 , with its unit 
vector basis {eo,ei}, let 

H m = H 1 ®---®Hi, 



equipped with the tensor Hilbert space structure. Denote 

Z[0,N) :={Q,...,N-1} 

for JV G N (we write N = {1, 2, . . . } and N = N U {0}). Let C m = {\i) : i G 
Z[0, 2 m )} be the canonical basis of H m , where \i) stands for ej ® ■ ■ ■ ® e Jm _ 1 
with i = Y1T=0 jk^ m ~ 1 ~ k · Let U{H m ) denote the set of unitary operators 
on H m . 

A quantum query on F is given by a tuple 

Q = (m,m',m",Z,T,P), 

where m,m' , m" G N, rrí + m" < m,Z Ç Z[0, 2 m ) is a nonempty subset, 
and 

r : Z -> A 
13 : K ^Z[0,2 m ") 

are arbitrary mappings. Let m(Q) := m denote the number of qubits of Q. 

Given a query Q, we define for each / G F the unitary operator Qf G 
U{H m ) by setting for |í) |x) |y) G 



Qf\i) \x) \y) 



\i)\x®(3{f{r(i))))\y) ü i£Z 
\i) \x) \y) otherwise, 



where © means addition modulo 2 m . 

A quantum algorithm on F with no measurement is a tuple 

A=(Q,(U j )] =0 ), 

where Q is a quantum query on F, n G No and Uj G U(H m ) (J = 0, . . . ,n) 
with m = m(Q). Given / G F, we define ^4/ G U{H m ) as 

A f = U n Q f U n - 1 ...U 1 Q f U . 
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We denote by n q (A) := n the number of queries and by m(A) = m = m(Q) 
the number of qubits of A. Let (Af(x, y)) x ,yez[0,2 m ) be the matrix of the 
transformation Af in the canonical basis C m . 

A quantum algorithm from F to G with k measurements is a tuple 

A = ((A l )1- \(b l )f- \ ( p), 

where k G N, A\ (l = 0, . . . , k — 1) are quantum algorithms on F with no 
measurements, 

b £Z[0,2 m »), 

2-1 

òj : JJZ[0,2 m< ) ->Z[0,2 m ') (1 <l< k-1), 

i=0 

where m\ := m(Ai), and 

fe-i 

¥>: JJZ[0,2 m; ) -»G. 
z=o 

The output of ^4 at input / £ F will be a probability measure -A(/) on G, 
defined as follows: First put 

PA./^o,---,^-!) = \A ,f(xo,bo)\ 2 \A l j(x 1 ,b 1 (x ))\ 2 ... 

• • • |4fc-l,/(x fe _iA_i(x , • • • ,«fe-2))| 2 - 

Then define A(f) by setting for any subset C Ç G 

J)(io,-A-i)€C 

Let := X^o 1 n <?(A) denote the number of queries used by A. For 

more details and background see [11]. Below we use the term 'quantum 
algorithm', meaning a quantum algorithm with measurement(s). Note that 
a quantum query on F (respectively, a quantum algorithm from F to G) can 
also be considered as a quantum query on ,F(A, K) (respectively, a quantum 
algorithm from J-(A, K) to G), and vice versa. 

The above definition simplifies essentially for an algorithm with one mea- 
surement, which is given by 

A=(Ao,bo,<p), A = (Q,(Uj)] =0 ). 
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The quantum computation is carried out on m := m(Q) qubits. For / € F 
the algorithm starts in the state |&o) and produces 

= U n QfU n - 1 ...U 1 Q f U \b ). 

Let 

2 m -l 
i=0 

Then the output takes the value ip(i) G G with probability |ai,/| 2 - As shown 
in [11], Lemma 1, for each algorithm A with k measurements there is an 
algorithm A with one measurement such that A(f) = A(f) for all / G F 
and A uses just twice the number of queries of A. 

For 9 > and a quantum algorithm A we define the (probabilistic) error 
at / G F as follows. Let ( be a random variable with distribution ^4(/). 
Then 

e(S,A,f,0) = inf {e > | P{||S(/) - Cll 
(observe that this infimum is always attained). Let 

e(S, A, F, 6) = sup e(S, A, f, 9) 
feF 

(this quantity can take the value +oo). Furthermore, we set 
e*(S,F,9) 

= inf{e(5, A, F, 9) \ A is any quantum algorithm with n q (A) < n}. 
We denote 

e(S,A,f )=e(S,A, /, 1/4) 

and similarly, 

e(S, A, F) = e(S, A, F, 1/4), e«(S, F) = e^S, F, 1/4). 

The quantity en(S,F) is the n-th minimal query error, that is, the smallest 
error which can be reached using at most n queries. Note that it essentially 
suffices to study en(S, F) instead of en{S,F,9), since with 0{y) repetitions, 
the error probability can be reduced to 2~ u (see Lemmas 3, 4 and Corollary 
1 of [14]). 

The quantum query complexity is defined for e > by 
comp^S, F) = 

min{n,j(^4) | A is any quantum algorithm with e(S,A,F) < e} 
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(we put compe(S, F) = +00 if there is no such algorithm). It is easily 
checked that these functions are inverse to each other in the following sense: 
For all n G No and e > 0, €n(S, F) < e if and only if compej (S, F) < n for 
all £1 > e. Hence it sufnces to determine one of them. We shall principally 
choose the first one. 

Note that the definition of a numerical problem we presented here cor- 
responds to that used in [17, 18] for the classical settings, and is slightly 
more general than the one in previous papers on the quantum setting [11, 
12, 14, 15]. There F was always a set of functions on some set D. We get 
back to this setting by considering, as done above, each / as a function on 
A and defining D = A. (Such an approach has already been outlined at the 
end of [11].) The mapping that sends / G F to the corresponding function 
(/(A))agA needs not to be one-to-one, in general. Nevertheless, all general 
results of [11, 12, 14, 15] carry over in an obvious way, with literally identical 
proofs. 

3 Some general results 

Let V = (F, G, S, K, A) be another numerical problem. Suppose we have 
an algorithm for problem V, and we want to construct one for problem V . 
Furthermore, for each input / G F of problem V we can produce an input 
R(f) for problem V such that S(f) = * o S o R(f) with a certain mapping 
^ : G — > G. Finally, each information about R(f) can be obtained from k 
suitable informations about /. Then we say that problem V reduces to V . 
Let us specify the assumptions. 

Let R : F — > F be a mapping such that there exist a k 6 N, mappings 
rjj : A — > A (j = 0, . . . , k — 1) and g : A x K K —> K with 

(R(f))(X) = q(X /Mà)), . . . , f(Vn-lW)) (i) 

for all / G F and A € A. Furthermore, let ^ : G —> G be a Lipschitz 
mapping and assume that 

5 = ^0^0^. (2) 

Note that (1) defines also a mapping 

R:F(A,K) -^T(A,K) 

(we use the same notation R), where ,F(A, K) stands for the set of all map- 
pings from A to K. 
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Lemma 1. Let Fq Ç F be any nonempty subset. Suppose that for each 
5 > and each finíte subset Aq Ç A í/iere are mappíngs 

6:K^K, G:F ^F 

such that 9(K) is a finíte set, 

(6(/))(A) = 0(/(A)) (/GF , AgAq), (3) 

and 

sup||S(/)-S(e(/))||<<5. (4) 

T/ten for all n G No, 

eLn(S,*b)< ||*||Li P eS(5,F). (5) 

Proof. Let í > 0, n G No and let vi be any quantum algorithm from F to G 
with ra ç (^4) < n and 

e(5,A,F)<eS(5,F)+(y. 

Let 

^4 = ((4)?- \ (^)f=o^) 5 4 = (Qí. (^j)"Io). 
= (mi,fh' h rh'í,Zi,Ti,Pi), 

where Ç Z[0,2™í) and 

tíiZí-A, À:^^Z[0,2™")- 

Denote 

Ao = {tí(í) : ieZ h l = 0,...,k-l} 

and 

Ao = {VjW ■ À G A }. 

Now let 6 and be according to the assumptions, and choose m* so that 
< 2 m *. It is easily checked that one can find 

P : K -^Z[0,2 m *) 

and 

7 : Z[0,2 m *) -► K 

such that 7 o (3 = 6. Define 

g:Ax Z[0,2 m *) K -► K 
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for Ae A, í ,...,í«_i eZ[0,2 m *) K by 

g(\, t ,..., í„_i) = £»(À, 7(t ), • • • , 7(*«-i)) 

and 

R:F{A,K) -*r(A,K) 

by 

R(f) = R(9of). 
Then for / G .F(A, AT) and A € A, 

(£(/))(A) = (fl(ío/))(À) = e (À,8o/( w (À)),.Jo/fe. 1 (À))) 
= ^(A,/3o/(r ?0 (A)),...,/?o/(r ?/t _ 1 (A))). 

Thus, the mapping R is of the form needed to apply Corollary 1 of [12]. 
Accordingly, considering A as a quantum algorithm from ^"(A, K) to G, we 
can find a quantum algorithm A from T(A,K) to G with n g (A) = 2Kn g (A) 
and 

A(/) = A(fí(/)) (fef(A,K)). 
For A G Ao we have r?j(A) G Ao, and therefore, by assumption (3), for f E Fq, 

(R(f))(X) = g(\, 6 o /(%(A)), ...,0o /(r/ K _i(À)))^ 

= ^(A,(e(/))(r ? o(A)),...,(6(/))(r ?(t _ 1 (A))) 

= (i?(6(/)))(A). 

This implies 

Ql,R(f) = Ql,R(<3>(f)) ( l = 0, . . . , , fe — 1), (6) 

and consequently 

= !(«(/)) = 2(i?(e(/))) (/ g F ). (7) 

Now fix / G Fo and let £ be a random variable with distribution We 
have, by assumption (4) 

IW)-*(0II < \\s(f)-s(e(f))\\ + \\s(G(f))-*(0\\ 

< \\S(@(f))-*(0\\+5. (8) 

Furthermore, by (2), 

||5(6(/))-*(C)|| = ||*oSofí(0(/))-*(C)|| 

< ||*|| Lip ||5oi?(e(/))-ci|. (9) 
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Since ©(/) e F,jwe have R(@(f)) e F. Moreover, by (7), the distribution 
of Ç is equal to A(R(@(f))). Therefore we get with probability at least 3/4, 

\\SoR(e(f))-Ç\\<e(S,À,F), 

and hence, by (9), 

||5(G(/))-*(C)|| < \\nu P e(S,I,F) 

< \\nu P K(S,F) + S). 

^(C) is a random variable with distribution ^f(A)(f) - the output of the 
quantum algorithm ^(A) from F to G (compare Lemma 2 of [11] and the 
definition before it), an algorithm with not more than 2/cn queries. This 
implies (5). □ 

We need some further notation. For a linear space X we denote by X* 
the algebraic dual, that is, the space of all linear (not necessarily continu- 
ous) functionals on X, and by X* the dual space, which is the space of all 
continuous linear functionals on X. Given a subset Fq of a normed space X 
and 5 > 0, we denote by Fq the closed 5-neighbourhood of Fq, that is, the 
set 

F$ = U xeFo B(x,S), 

with B(x, d) being the closed ball of radius 5 around x. The unit ball B(0, 1) 
of X is denoted by F>x- 

Lemma 2. Let K = 1K, let F be a bounded subset of a normed space X , 
and let / Fq Ç F. Assume that either 

(i) there is a So > such that F^° Ç F or 

(ii) F is a non-zero múltiple of the unit ball of X. 
Furthermore, let Ao C X# be a finite, linearly independent set with 

sup feFo \f(X)\ < oo (Xe A ). 

Then for each 5 > there are mappíngs 

6 : X -► X, 9 : F -► F 

suc/í í/ioí 6(K) is a finite set, 

(e(f))(X) = 6(f(X)) (f€F , AgAq), (10) 

anc? 

su P ||/-e(/)||<í. (ii) 

/eF 
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Proof. We can assume S < 5q < 1. The linear independence of Ao implies 
that for each A G A there is a g\ G X with g\{\) = 1 and g\(n) = for 
fi e A \ {A}. Define 



Mi = max|| 5A ||, M 2 = sup |/(A)|, M 3 = sup 
AeA AeA ,/GF /eF 

íi = <5/(M 3 + 1), 

and choose any #o : K — > K such that 0o(K) is finite and 

\a - 9 (a)\ < Mf^Aol -1 *! min(M 3 , 1) (|o| < M 2 ). 

Now we define 6 : K — > K by setting for oGl, 

9{a) = e ((l-5 1 )a), 

and 6 : F -> X by 

G(/) = (1 - <5i)/ - (í 1 - <W( A ) " ^o((l " 5i)/(A)) 5A . 
AeA 

Then for / G F , \i G A , 
(©(/))(M) 

= (1 - <5i)/(/x) - ^ ((1 - <5i)/(A) - 0„((1 - <5i)/(A)) 5a (m) 
AeA 

= (9 ((1-Íi)/(m)) = Wa*)), 
which verifies (10). Moreover, we have, by (12) and (14), 



2((l-íi)/(A)-e ((l-íi)/(A))^ 

AeAo 



(12) 
(13) 

(14) 



(15) 



< 



£ l(l-íi)/(A)-flb((l-íi)/(A))|í/ A 

AeA 



< min(M 3 ,l)(5i. (16) 



Hence, by (15), (16), and (13) 

11/ - 6(/)|| < b x \\f\\ + <5i < (M 3 + l)íi = 5 < So, 

which proves (11). Furthermore, it shows that in case of condition (i), 
0(/) G F for all / G F . If condition (ii) is fulfilled, that is, F = a B x for 
some ao > 0, we argue as follows: 



||e(/)|| < ||(1 - + M 30l < (1 - ó^ao + íia = a , 

thus, again, ©(/) G F for all / G F . 



□ 
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Proposition 1. Let K = K. Assume that S,S,R,^> are as above (1), 
(2), that F is a bounded subset of a normed space X, and A is a línearly 
independent subset o/X*. Let Fq be a nonempty subset of F and assume 
that eíther 

(i) Fq° Ç F for some 5 > 0, or 

(ii) F is a non-zero múltiple of the unit ball of X. 

Furthermore suppose supj g ^ |/(A)| < oo for each A G A and S is uníformly 
continuous on F . Then for all n G No, 

e1 Kn (S,F )< \\*\\u v <%(S,F). 

Proof. This is a direct consequence of Lemmas 1, 2, and the uniform conti- 
nuity of S. □ 

In previous papers on quantum complexity [11, 12, 15] the analysis of 
reductions was somewhat cumbersome, since a certain discretization had to 
be applied in each particular case. Proposition 1 simplifies the analysis and 
will be used for a number of reductions, in particular in sections 5 and 6. 

Next we recali additivity properties of the quantum minimal error, see 
[12], Corollary 2. 

Proposition 2. Let p G N and let Si : F — > G (l = 1, . . . ,p) be mappings 
such that S(f) = Y^d=i^l(f) (/ G F)- Let ui,...,u p <E N be numbers 
satisfying 

p 



i=i 



Then for all n\, . . . , n p € No 

=1 



Given a subset B Ç X of a normed space X, we denote by 'té'(B) the set 
of all precompact subsets of B. A set H C X* is called linear ly independent 
over a non-empty set B Ç X, if the restrictions of elements of H to span(S) 
form a linearly independent subset of (span(S))#. 

Finally we state multiplicativity properties of the minimal quantum er- 
ror. 

Proposition 3. Let K = K. Assume that F is a subset of a normed space 
Y , that A is a linearly independent subset of and supy g ^ |/(A)| < oo 
for each A G A. Leí J : F — > V òe í/ie embedding map, let T : Y — > G 
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be a bounded linear operator and assume that S = TJ. Furthermore, let 
v\ , V2 £ N be any numbers with 

e -"i/8 + e -»2/8 _ e -(vi+v 2 )/8 < 1 j Am 

Along with V = (F, G, S, K, A) we consider the problems (F, Y, J, K, A) and 
(B Y ,G,T,K,A). Thenfor alln 1 ,n 2 G N , 

4 ini+ 2u 2 n 2 (S,F) < 4e« 1 (J,F)e« a (T,B Y ). (18) 

//, moreover, F is a precompact subset of Y and A is linearly independent 
over F, then 

el ini+2u2n2 (S,F)<4el(J,F) sup e^(T,S). (19) 

The first part, relation (18), was proved in [14], Proposition 1 and Corol- 
lary 3. As already mentioned at the end of the previous section, this result 
was formulated for a slightly less general type of numerical problem, but the 
proof of (18) is literally the same as that of Proposition 1 and Corollary 3 
in [14]. 

The specific form of multiplicativity stated in (19) (a single en{S, F) 
is replaced by a supremum over a family of subsets of F) will be needed 
in sections 5 and 6. For further explanation we refer to the remark after 
Proposition 6. 

Proof of Proposition 3. It remains to prové the second part, relation (19). 
We derive it from the first part, (18). Denote 

e^(J,F) = a (20) 

and fix any 5 > 0. Let A = ((^^=0*' (Pi)i=o^) ' 3e a <l uan tum algorithm 
from F to Y with n q (A) < n\ and 

e( J, A,F) < a + 5. (21) 

Let Ç be a random variable with distribution A(f). Observe that, by defi- 
nition (see section 2), £ takes vàlues in the finite set 

y = ^nz[o,2"")j CY. 

Define 6 o C Y to be the closed, absolutely convex hull of F U Yq, and put 

£ = B Y n -^— x £o- 

a + ò 
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Since FL)Yq is precompact in Y, so are £q and £. Moreover, there is a 7 > 
such that F UY Ç 7Fy. Hence, 

F U y Ç max (7, £. (22) 



For any / G F we have 

/ - C e 2£ 
and by (21), with probability at least 3/4, 

/-Ce (<7 + í)fly. 

Consequently, with probability at least 3/4, 

/-(6(<r + 6)8. (23) 

Since £ is a closed, absolutely convex, and bounded subset of Y, it defines 
a norm || . ||g on F = span(£) as follows 

|| y || £ = inf{0 > : y G 05} (y £ E), 

and 5 is the unit ball of (F, || \\ £ ). By (22), F C E and F C F. Define 
Je ■ F E and ^ to be J and <£>, respectively, considered as mappings into 
E. Define A E = ((A,)^ 1 , (ò^t q^e)- Then A# is a quantum algorithm 
from F to E with íi g (A B ) < m. By (23) and (20), 

< ( J E , F) < e( J £ , A £ , F) < <j + í = e\ ( J, F) + 5. 

Note that, since F Ç F, A is linearly independent over F. Furthermore, 
since Be = £ is bounded in Y, the restriction of T to F is a bounded linear 
operator from F to G. Applying now the first part of Proposition 3, we get 

<n 1+ 2„ 2 n 2 (^) < 4e* 1 {J E ,F)e« a {T,B E ) < 4{e« 1 {J,F)+6)e% t {T,S), 

which gives the desired result, since 5 > was arbitrary. □ 

Let e^ ct (S', F) denote the (classical) n-th minimal deterministic error, 
that is, the minimal error among all deterministic, adaptive algorithms using 
at most n informations (see, e.g., [17], section 4). We want to apply relations 
(18) and (19) with ( J, F) replaced by e£f ( J, F). An estimate of e&(S, F) 
by e^ et (S', F) is not obvious, since classical deterministic algorithms can use 
information with vàlues in K directly, while quantum algorithms can use 
them only through a finite encoding. We therefore supply the following 
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Lemma 3. Let K = K, assume that F is the unit ball of a normed space 
X, A Ç X* , and S is a bounded linear operator from X to G. Then for all 
n £ N 

e*(S,F)<2e??(S,F). (24) 

Proof. Let 5 > 0. It is well-known (see [28], Theorem 5.2.1 and Corollary 
5.2.1) that there is a nonadaptive deterministic algorithm A using at most 
n informations such that 

svp\\S{f)-À(f)\\<2e% st {S,F)+6. (25) 
feF 

Such an A has the following form: There are Ao, • • • , A n _i G A and a mapping 
ip : K n -> G such that 

A(/) = ^(/(A ),...,/(A n _i)) (feF). (26) 

Without loss of generality we may assume that Ao = {Ao, . . . , A n _i} C X* is 
a linear ly independent set (if not, we pass to an independent subset and omit 
the rest by suitably modifying (p) . Let 6 and 6 be the mappings which result 
from the application of Lemma 2, case (ii). Put m' = 1, choose m" G N such 
that |6*(K)| < 2 m and let m = m! + m" . We represent (as done before, in 
the proof of Lemma 1) 

with (3 :K ^ Z[0,2 m ") and 7 : Z[0, 2 m ") -> K. Furthermore, we identify 
Z[0, 2 m ) with {0, 1} x Z[0, 2 m "). Define ^ : (Z[0, 2™))™ G for 

(6i, Oi) G {0, 1} x Z[0, 2 m ") = Z[0, 2 m ) (i = 0, . . . , n - 1) 

by setting 

<p((& ,a ), • • • , (6 n _i,a n _i)) = 9^(7(00), • • • ,7(on-i))- 
Now we construct a quantum algorithm A with n measurements. We let 

A=((A)r= 1 ,w=o 1 ^)> 

with 60 = 0, b\ = (1 < l < n — 1), and (p as above. Each A\ is of the form 

M = (Qí,(^)j=o,i) 
with Z7o = í/i = lH m the identity matrix, 

Qi = (m,m',m",Zi,T,/3), 
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with m, m', m", [3 as defined above, Z\ = {0} and r/(0) = Aj (Z = 0, . . . , n— 1). 
This simply means that A is an algorithm which queries the function / 
in the appropriate n points, with the needed precision, and measures the 
result after each query. Finally <p is applied. Let / G F and let Ç have 
distribution A(f). Since the computation remains on the classical states, 
the measurements give the result 

(0,/3(/(A ))),..·,(0,/3(/(A n _ 1 ))) 

with probability 1. Hence 

C = £(/W(Ao),...,/W(A„-i)) 

= v 3 (7°/5 /(Ao),···,7 /3°/(An-i)) 
= ^o/(A ),...^o/(A n _ 1 ))=i(6(/)) 

with probability 1, consequently, by conclusion (11) of Lemma 2 and by (25) 

IW)-cil < \\s(f)-s(@(f))\\ + \\s(@(f))-À(e(f))\\ 

< \\S\\S + 2et ct (S,F)+5, 
and (24) follows. □ 

From Proposition 3 and Lemma 3 we immediately conclude 

Proposition 4. Let X and Y be normed linear spaces such that X is a 
linear subspace of Y and the embedding J : X — > Y is continuous. Assume 
that K = K, F = Bx, A is a linearly independent subset ofY*, T :Y — > G 
is a bounded linear operator, and S = TJ. Then for all m, n2 G No, 

^n 1+ 2u 2 n 2 (S,F) < 8e«f(J,F)e« 2 (T,B Y ), (27) 

where u\, vi are any numbers satisfying (17). If, furthermore, J is a compact 
operator and A is linearly independent over X , then 

^ 1 „ 1+2 , 2 n 2 (^^)<8e^ t (J, J F) sup e^(T,S). (28) 



4 Weighted mean computation and integration 

Let Li, respectively, L^, be the space of all functions / : Z[0, N) —>■ K, 
equipped with the norm 

1 N - 1 

11/11^ = ^21/(01. 

i=0 
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respectively, 



t n = max \f(i)\. 



We use the notation L±(K) and L^(K) if the underlying field has to be 
emphasized. Let g € L^. Define the weighted mean operator 5jv, s : — ► 
K by 

1 iV " 1 
i=0 

We write S^v for Sn 9 with g = 1. We consider the weighted summation 
problem p = (B L N,k,S N , g ,K,A) with A = {Si : < i < N} and $(/) = 
/(i). Throughout this paper we often use the same symbol c,c\,... for 
possibly different positive constants (also when they appear in a sequence 
of relations). 

Proposition 5. There is a constant c > such that for all n, N G N and 

4(SV, 9 ,£ L ^) < cn ^{çWln. 

Proof. First we consider the case K = R. If g = 0, the statement is trivial. 
We may assume without loss of gener ality that g > 0, otherwise we split g 
into its positive and negative part and apply Proposition 2. Moreover, by 
scaling the problem appropriately, we can assume 

Ml» = !■ ( 29 ) 

Now we reduce the problem Sn,ç to the known case Sm for some M. Define 
h,g€L?by 

h(i) = [ng(i)\, g(i) = n~ 1 h(i) (i = 0, . . . , TV - 1). (30) 

We have \g(i) — g(i)\ < n -1 , therefore 

sup \S N: gf - S N> gf \ < n -1 . 
feB jv 

By Lemma 6 of [11], 

<e^(5 JV , 5 ,B i ^) + n- 1 . (31) 
Now set mo = and for 1 < i < N 

i-i 

m = EM0) 

Z=0 
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and denote itin = M. The case M = is trivial, since this implies g = 0, 
thus the result follows directly from (31). Hence we assume M > 1. Observe 
that (29) and (30) imply 

M < nN. (32) 

Define 

7] : Z[0,M) -► Z[0,iV) 

by = i, where i is the unique integer satisfying m; < j < mj+i. Let the 
reduction mapping i? : — > be given by 

WW) = f(.vU)) ü' = o,...,m-i). 

Clearly, i? is of the form (1), with k = 1, therefore, by Proposition 1 we 
have 

e* n {S M R,B L ^) <el{S M ,B L M). (33) 

Moreover 

M-l AT-l 

s M R{f) = ü E /(i/ü)) = m E 

j=0 i=0 

= mE^)/(í) = ¥ %/. 

This together with (32) and (33) implies 

$ n (S N ,g,B LS) ) = e% n (J^S M R,B L N^ =-^eï n (S M R,B L N) 
< e%(S M ,B L M) < cn -1 , 

the latter relation being a consequence of [3] (see also [11], Theorem 1, for 
the form stated here). Combining this with (31) and scaling the index gives 
the desired result. 

Now we formally derive the complex case from the real case. Let g € 
Lf (C) and let gi ,g 2 € (R) be defined by 

gU)=9iü) + mU) (j = o,...,N-i, i = V=ï). (34) 

Clearly, 

||5a||Lf( K ) < llflLf (C) (« = 1 > 2 )· ( 35 ) 
We shall express V = (B l n , C, Sjv, ff , ^ A) by the help of 
P=(B L£ ,R,,S' JVií , a ,K,A) (a = l,2). 



17 



Define R ± ,R 2 : i£(C) - L£(R) for / G L£(C) by 

(R 1 f)(j)=Re(f(j)), (i? 2 /)(j)=Im(/(i)) (i = 0, . . . , N - 1). (36) 

Clearly, R\,R2 are of the form (1) and map B l n^ to B l nç^. Define 
J al3 :R^C (a,0€ {1,2}) by 

Jua = —J22C1 = a, J\2ü = J21C1 = ta (a G M). 

Then we have, by (34) and (36), 

2 

S N,gf = ^2 J aP S N,g a R pf- 

Let v be the smallest natural number with e~ u / 8 < 1/16. By Proposition 2 

2 

e 2un (S N;g ,B L N( C · ) ) < 2 ^ e2n(Jaf3S N ,g a R / 3,B L N( C - ) ). (37) 

Q,/3=l 

Moreover, by Proposition 1, 

e2n{JapS N>ga Rp,B L N(c)) < en(S Ní g a ,B L N( R )) (a,/3 = 1,2). (38) 
Using the result for the real case and (35), (37), and (38), we get 

e2vn(SN,g,B L N( C )) < crT 1 \\g\\ l n (c) , 
and a scaling of the index concludes the proof. 

□ 

Now we pass to the case of weighted integration. Let Q Ç M. d be a 
closed, bounded set of positive Lebesgue measure. L±(Q) denotes the space 
of Lebesgue integrable functions on Q with vàlues in K, equipped with the 
norm 

II/IImq)= / 

JQ 

and Loo(Q) the space of all K-valued measurable and essentially bounded 
with respect to the Lebesgue measure functions on Q, endowed with the 
norm 

H/lloo = esssup x6Q |/(x)|. 
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Let g G L\(Q). Define Iq j9 : L 00 (Q) — ► K, the integration operator with 
weight g, by 

1 qJ= \ g{x)f(x)dx. 
JQ 

C(Q) denotes the space of continuous functions on Q, equipped with the 
supremum norm. A set £ of continuous functions on Q is called uniformly 
equicontinuous, if for each e > there is a 5 > such that for x,y G Q, 
\x — y\ < 6 implies \f(x) — f(y)\ < e for all f E F. By the Arzelà-Ascoli 
theorem, bounded, uniformly equicontinuous sets coincide with precompact 
subsets of C(Q). We consider the problem V = (Bç^, K, Iq :9 , K, A) with 
A = : x £ Q}, where í x (/) = /(x) for / G C*(Q). 

Proposition 6. There is a constant c > such that for each closed, bounded 
set Q CR d of positive Lebesgue measure, for all g G L\(Q) and n G N 

f sup (/q,í„ f) < cn _1 ||5|| Ll ( Q) . 

£e c <f(B C( Q)) 

Remark. It is well-known and easily checked by using importance sampling 
with density function IsÍ/HpIIl^q) that in the classical randomized setting 
we have 

e^{I Q , g ,B c{Q) )<cn-^\\g\\ Ll{Q) , 

where e™ 11 is the n-th minimal classical randomized error (see, e.g., [18], 
section 3). Proposition 6 is the quantum analogue of this result. Let us 
comment on the reasons for taking the supremum over £ G ^{Bçiq))- In 
contrast to the classical randomized setting, no non-trivial convergence rate 
holds for en(lQ t9 , B c ^), in general. This is easily checked based on the fact 
that a quantum query involves, by definition, the vàlues of functions from 
Bc{Q) in a finite set of points of Q only. For situations like this a natural way 
of formulating quantum counterparts of results of the classical randomized 
setting was already observed in section 5 of [11]: If we restrict our analysis 
to uniformly equicontinuous subsets £ of the respective unit ball, non-trivial 
decay rates can be shown in such a way that neither the exponent nor the 
constants involved in these estimates depend on £ (though the number of 
qubits in the respective quantum algorithms does, but this is irrelevant for 

f&(lQ,g,S)). 

Proof of Proposition 6. Fix Q C R d , £ G C €(B C(Q) ), and n G N. Let Q* be 
a cube with Q Ç Q*. For k G N let 

2 dfc -l 

Q*= U Q*i 

i=0 
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be the partition of Q* into 2 dk congruent cubes of disjoint interior. Let 
Qi = Q* H Q- Without loss of generality we assume them ordered in such a 
way that n(Qi) > iff i < N, where n is the Lebesgue measure and ./V is an 
appropriate number 1 < N < 2 dk . Then 

jv-i / N-i \ 

[J Q t Ç Q and n Q \ (J Qi = 0. 

i=0 V i=0 / 

Let Xi be any point in Qi and let be the operator of piecewise constant 
interpolation with respect to the partition (Qi)^S^ 1 in the points (xj)^ 1 . 
By the uniform equicontinuity of £, there is a k such that 

\\f-I\f\\L oaí Q)<n- 1 (39) 

for all f £ £. It follows that 

sup\I Q:9 f - I Q:g (P k f)\ < n^WgWL^qy (40) 
fes 

We define R : B C(Q) - by 

= /(xí) (i = 0,...,JV-l). 

Then i? is of the form (1) with A = {Si : < i < AT} and maps #c(Q) to 
-B l at . Furthermore, define /i G by /i(i) = N Jq g(y)dy. Then 



and 

N-l 



lQ, 9 (Pkf) = E /(^) / 9(t)dt 
1 N ~ 1 

= ^EMi)W))(i) = 5jv >h ol2(/). (41) 

j=0 

Lemma 6 of [11] together with relations (40) and (41) imply 

e*{I Q , g ,£) < n^Wgh^ + e*(S N , h o i2,f ). (42) 
By Propositions 1 and 5, 

4n( S N,h°R,£) < e n( S N,h,B L N) < en' 1 \\ h\\ L N < en' 1 \\g \\ Ll (Q) , 
which together with (42) accomplishes the proof. □ 
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5 Quantum approximation of weakly singular in- 
tegral operators 

Let 1 < di < d and let Q\ be the closure of an open bounded set in R dl . We 
identify Qi with a subset of R d by identifying R dl with M dl x {O^ 1 )}. Let 
Q2 be a bounded Lebesgue measurable subset of R d of positive Lebesgue 
measure and define diag(Qi, Q2) := {(x,x) : x G Q\ n C^}- We introduce 
the following class of kernels (see also [17], where integral operators with 
such kernels are analyzed). 

For s G N and (ieR with -d < a < +00 we denote by C s ' a (Qi, Q 2 ) the 
set of all Lebesgue measurable functions k : Q\ x Q2 \ diag(Qi, Q2) — ► K 
with the following properties: There is a constant c > such that for all 
V G Q2 

1. k(x,y) is s-times continuously differentiable with respect to x on 

Qi \ {y}' where Q\ denotes the interior of Qi, considered as a subset 
of R d \ 

2. for all multiindices a G N dl with < |a| = a± + ■ ■ ■ + < s the 
a-th partial derivative D"k(x,y) of k with respect to the x-variables 
satisfies the estimate 

for all x £ Qi\ {y}, and 

3. for all a G N dl with < \a\ < s the functions D"k(x,y) have contin- 
uous extensions to Q\ \ {y}. 

We want to extend the definition to the case d\ = 0. Here we let Q\ = 
{0} C R d and define C s ' a (Q 1 ,Q 2 ) to be the set of all functions fc(0, y) which 
are Lebesgue measurable in y and satisfy 

l*(0,y)|<c{ + 1 !| IÍ° (yeQ 2 \{0}) (44) 

with a certain c > 0. Note that for di = the space C s,cr (Qi, Q2) does not 
depend on s. 

For k G C s ' CT (Qi,(52) let ||fe||c«,<7 be the smallest c > satisfying (43) or 
(44), respectively. It is easily checked that || . \\c«," is a norm on C s,<T (Qi, Q 2 )- 
For G C s,a (Qi, Q2) we let be the integral operator 

(T k f)(x)= í k(x,y)f(y)dy (x G Qi) 
JQ2 
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acting from C{Q2) to Loo(Qi) (to K, if di = 0). We shall also consider 
Tfc as acting in various other function spaces, which will then be mentioned 
explicitly. It is easily checked that maps C(Q 2 ) into C(Qi). Finally, 
denote 

C°°' a (Q 1 ,Q 2 ) :=f)C s '°(Q 1 ,Q 2 ). 

seN 

We start the analysis with the case of Q\ = [0, l] dl , where < d\ < d, and 
Q2 being a closed subset of [0, \} d of positive Lebesgue measure. We study 
the minimal quantum error of approximating T^f with k G C s,CT (Qi,Q2) a 
fixed kernel, thus we consider F = B C (q 2 ^, G = L^Qi), S = T^, K = K, 
and A = {5 X : x G Q2} (for d\ = the space L^Qi) is replaced by K). 

To state the following proposition, define (3{cr) (this parameter will de- 
scribe the power of the logarithmic term) as 

if min(s, d + a, d) > d\ 
4 if min(s, d + a,d) = d\ 

= { min(s,d+a) -r ■ / , . x . , , , , . ( 45 ) 

— ^ '- if mm(s, d + <j) < di and s 7^ d + <r 

min ^' i ci+cr · > +1 if min(s, d + a) < ííi and s = d + a. 

Note that, since d\ < (i, we have min(s, d + a) < d\ iff min(s, d + o~,d) < 
d\, so (45) covers all possible cases. The following is the quantum version 
of Proposition 1 of [17]. For the appearance of the supremum over £ G 
^{Bc^q^) we refer to the remark after Proposition 6. In the case di = 
we intèrpret ^ = = +00. 

Proposition 7. Zeí < di < d, s G N ; íj G E ; -d < a < +00, Qi = 
[0,l] dl . T/ten í/iere ís a constant c > siíc/i í/iaí /or any closed subset 
Q2 Ç [0, l] d of positive Lebesgue measure, and for all k G C s '°"((5i, Q2) 
n G N with n>2, 

sup e£(T fc ,£:) < C rT min (t'^f - 1 ) (logn)^)||fc|| c .,. ( Q 1)Qa) , (46) 
where (3(a) is as defined in (45). 

Proof. In view of Lemma 6 (ii) of [11] it sufnces to prové the statement for 
k with 

\\k\y.°( Qu Q2) = !• (47) 

In the case di = we have 



T k f= í k(0,y)f(y)dy, 

JQo. 
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where by (47) and (44) 

/ \k(0,y)\dy<c 
JQi 

(the constants in this proof depend only on d, d\, s, a), and the result follows 
directly form Proposition 6. 

Now we assume d\ > 1. First we recali some notation from [17]. For 
1 = 0,1,... let 

Qi = (j Qi,u (48) 
1=1 

be the partition of Q\ into ni = 2 dl1 closed subcubes of sidelength 2~ l and 
mutually disjoint interior. Let T[ be the equidistant mesh on Q\ with mesh- 
size 2"'(max(s - 1, l)) -1 , Tu = T/ n Qi t u and Tu = T i+1 n Qi t u. Let E u be 
the subspace of C(Qi j u) consisting of all multivariate polynomials on Q 1 j i 
of degree at most max(s — 1, 1) in each variable. Let Ei be the respective 
space of continuous piecewise polynomial functions on Q\, that is 

£?i = {/ € C(gi) : f\ Qhli eE H , i = i,...,m}. 

Furthermore, define En C C(Qi j u) by 

Eu = {/ € C(Q hU ) : f\Q hl+1J € E i+ld for all j with Q h i +hj C Q hU }, 

in other words, En C C(Qi t u) is the space of continuous piecewise polyno- 
mial functions with respect to the partition of Q\ t u into subcubes of side- 
length 2~( l+l \ Let Pu : /<x>(Fh) — ► En be the multivariate (tensor product) 
Lagrange interpolation of degree max(s — 1, 1) on Tu, define Pu : loo(ï"n) — > 
Eu by 

( p ii u )\Qi,i+i,3 = p i+i,j( u \r l+1 J 
for all j with C Qi^, and Pi : /oo(F/) -> E x by 

^«IQi,k = ^iMr K ) (i = l,...,ni) 

(continuity follows from the assumption that the degree is > 1). Thus, Pu 
and P\ are piecewise Lagrange interpolation operators. For / £ C(Qiu) 
or / £ C(Qi) we write P^/ instead of Pu(f\r u ), and similarly P^/ and 
Pif. We shall use the following well-known (see, e.g., [4]) properties: For all 
l € No and i = 1 , . . . , ni , 

\\Pu : loo(T H ) ^ C(Q 1M )\\ < c, (49) 
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furthermore, for / G C s (Qi t u), 

\\f ~ ^f\\c(Q 1M ) < c2- sl y\\ C s (QlM) , (50) 

and consequently, 

11(4 - Pii)f\\c(Q 1M ) < c2- s '||/|| Cs(QlH) . (51) 

Define the embedding operators Ju : C(Qi : u) — > L OQ (Qi) by setting for 
x G Qi, 

(Jlif)(x) | ^ ( ) otherwise. 

We identify C{Q\) with a subspace of Loo(Qi)> thus, the operators Pi can 
also be considered as acting into Loo(Qi). First we approximate T k f by 
PmTkf, where m > 1 will be fixed later, then P m T k f will be approximated 
by a quantum algorithm. It is readily checked that 

m— 1 n; 

PmTkf = P T k f+Y,Y, Jli( · pH - P ^ Tk f- ( 52 ) 

For / = 0, 1, . . . to— 1 and i = 1, . . . , n\ let be the center and q\ = \fd[2^ l ^ x 
the radius of Qi.zí- For £ > let B(x,g) denote the closed ci-dimensional 
ball of radius g around x G M. d . We represent 

(T fe /)(x) = / k(x,y)f(y)dy 
J B{x u ,2 ei )f\Q 2 

+ í k(x,y)f(y)dy, (53) 

jQ2\B(xu,2 ei ) 

and introduce kn G C(Qi : u x (Q2 \ B(xu,2qi))) by setting for y G Q2 \ 
B(xu, 2q{) 

k li (-,y) = (P li -P li )k(-,y). (54) 
Using that (P K - P K ) = (Pi - Pi) 2 ; we conclude 

Ju(Pu - Pu)T k f = MPii-Pu) í k(-,y)f(y)dy 

J B(x li ,2 Bl )nQ 2 

+Ju{Pu-P u ) í k li (-,y)f(y)dy. (55) 

J Q 2 \B(x lu 2 ei ) 

Next we introduce the following functions: For x G To define G L\(Q2) 
as 

k x {y) = k{x,y) (yeQ 2 ), 
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and for l = 0, . . . , m - 1, i = 1, . . . , n h x G T u define gu x , h hx G L 1 {Q 2 ) for 
y G Q 2 by 



9ux(y) 
hux(y) 



k(x,y) if y e B(x H ,2g l )nQ 2 
otherwise, 

hi{x,y) if y e Q 2 \B(x u ,2qi) 
otherwise. 



Thcn 



P m Tkf = Po ((lQ 2 ,k x f)xer ) 

m—l ni 

+ ^2Y1 Jh ^ h - p ^ (vq^uJ + / 02,hHj)o : 6f ii ) • ( 56 ) 

1=0 i=l 

From (47) and (43) we have 

IIMl i( Q 2 )=/ \k(x,y)\dy<c (x G r ). (57) 
For x G f/j we deduce from (43) that 



HfitoxIUiNa) = / |A:(ar, j/) | dj/ 

J B(xn,2 ei )nQ2 

< / |fc(x,y)|efy 

•/B(x,3ej)nQ2 

< c/ (|x-y| a + |ln|x-y|| + l)(íy 
JB(x,3oi) 



B(x,3ei) 

< c{2-( d+a V + (l + l)2- dl ). (58) 
Furthermore, again by (43), we have for x G Qi,u and y £ Q 2 \ B(xu,2qi) 



\D^k{x,y)\ < cdx-y^-W + llnlx-yll + l) (59) 
< c(\x lt -yr^ + \\n\x tí -y\ | + 1), 

hence 

|| K ■ , y)\\c°{Q 1M ) < c {\x u - yr s + | ln \x u - y\ | + 1). (60) 
Using (51) and (54) we obtain for y G Q 2 \ B(xu,2qi) 

\\ki( ■,y)\\c( Ql , u ) < c2- sl (\x u - y\°- s + | ln \x u - y\\ + 1). (61) 



25 



We have 

. fl iï a — s > —d 

/ \xu-y\ a ~ s dy < c< l + l if a - s = -d 

Jq 2 \b(xh,2qi) j 2-(< 7 - s + d ) 1 if a — s < -d. 

Therefore, integrating (61), we get for x G tu 

ll^lk(Q 2 ) = / \ku(x,y)\dy <c(l + i)«o 2 -min( s , d+CT )^ (g2) 

jQ 2 \B(x H ,2 m ) 

where 

a ={ 1 Íf s = d + a (q 3 ) 
\ otherwise. 

Now we approximate the integrals in (56) by quantum algorithms. Let 
£ G ^(Bq^q^) (as already mentioned, the constants depenc! only on the 
parameters d,d\,s,a, and in particular not on £). Using Proposition 6 
together with (57), (58), and (62), we obtain the following relations 

e q Nm VQ2,k x ,£)<c/NW (xer ) (64) 
e%(lQ 2 , gilx ,£) < c(2"( rf+ ^ + (/ + l)2- dl )N- 1 {x G f„) (65) 
*%VQ>to.J) < C V + ir^-^^N- 1 (x G t H ), (66) 

where iV* ), iVj G N (/ = 0, . . . , m — 1) are arbitrary natural numbers which 
will be fixed later. Let ,v\ G N (l = 0, . . . , m — 1) be the smallest natural 
numbers satisfying 

| r (o)| e -,(°>/8 < 2 -3 (67) 

ni 

2^|f /í |e- iy '/ 8 < 2-C+ 4 ) (Z = 0,...,m-1). (68) 



i=i 

Consequently 

^<c(Z + l) (1 = 0, ...,m-\) (69) 

and 

m— 1 ni 

| r (0)| e -^)/8 + 2 \t H \e-^ 8 < 1/4. (70) 

1=0 i=l 

By Lemma 3 of [11] we can assert the existence of quantum algorithms 
A ( x 0) (x G To) with 

e(/Q 2 , fc „4°U,e- (ü) /8)< c /Ar(0), ^W) < *,(°)jv(°) (x g Tq) (71) 
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and of quantum algorithms and (x € Tu, i = 1 . . . nj, l = 0, . . . , m— 
1), such that 



e^^^S^.^ 78 ) < c(2-(*^)' + (/ + l)2^)iV r 1 



and 



e^^.^S.f ,e-^ 8 ) < c(Z + l^-minM+^-i 
for all x G IV We define the quantum algorithm 

m— 1 n; 

Jh(Ph-Ph) í ía\}]+a\T, 

xer u 



(72) 
(73) 

(74) 



/ \ m—í ni , n. 

* ^ p » (K'L J + E E *A - ((42 + < 75 > 



in the sense of the composition of quantum algorithms described in [11], 
relation (11). Let / € S and let C*°\C&\C,ÏÏ G r K , i = l...n,, / = 
0, . . . ,m — 1) be independent random variables with distribution A x °\f), 
A^2(/)' anc ^ respectively. From Lemma 2 of [11] it follows that the 

random variable 

(v m-1 ni , v 

^U + EE^-«.)((cS + cS) ie J (76) 

has distribution and that 

m— 1 n 

n q (A) = ^ ' " nA 

x6 r(°) z=o í=í xeTl 

By (56), (76), and (49), 
l|ïfc/-CII 



m— i n; 

= £ (4 0) ) + + n 9 (42)) (77) 



< \\T k f - P m T k f\\ + 

TO— 1 «I 



+ 



EE - p ià ( -/ - <£!) 



Z=0 i=l 



ser,) 



+ - -Pm) ( {lQ2,h lix f ~ Cul) 



xer H 



< \\T k f - P m T k f\\ + cmax I Q ^J - 



TO— 1 TLl 

+ c EE 

Z=0 i=l 



max 



In f - 



+ max 



1 Q2,hiixJ ^lix 



(78) 
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As established in [17], Lemma 3, the error of approximation by P m Tkf 
satisfies 

\\T k f - P m T k f\\ < c ( m <*o 2 -minM+a)m + m2 ~ dm ). (79) 

Furthermore, from (70), (71), (73), and (74), we conclude that with proba- 
bility at least 3/4 the following relations hold simultaneously: 

max\I Q2 ,kJ-Ú 0) \<c/NW (80) 

and for i = 1 . . . n/, / = 0, . . . , m — 1 

max \I Q2 , gu J - C«| < 0(2-^)' + (Z + l)2^)iV r 1 (81) 

and 

max |/ Q2Ai J - Cg| < c(l + i^-minM-H^-^ (g2) 

We get from (78-82) 
e(T fc ,A,f) 

< c ( m «0 2 -min( S ,d+ CT )m + m2 -'im) + c (jV(°)) _1 + 
m— 1 

c iVf 1 ((/ + l) ao 2" min (M+<^ + (/ + 1)2-*) . (83) 
By (69), (71), (72), and (77), the number of quantum queries of A satisfies 

m—í ni 



n q {A) < ^ (0) iV( )|r( )| + 2^^^iV / |f /í | 

1=0 i=l 

(m-l \ 
iV(°) + ^(/ + l)2 dlí iV í . (84) 

Let n G N with n > 2 be given. First we consider the case min(s, d + a, d) > 
di. Let r > be such that 

min(s, d + a,d) > di + r, 

and put 



m 



logn 
Ji + t 
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(log always means log 2 ), 

AT(°) = n , Ni 



n2 



-(di+r)Z 



(1 = 0,..., m-l). 



By (83) we get for the quantum error: 



e(T fc ,A,f) 

< c ( m °o 2 -min( S ,d+a)m + m j-rfm + n ~l) 
m— 1 

+C ^ n- 1 2( dl+T )' ((/ + l)«0 2 -min( S ,d+a)i + ^ + 1 ) 2 ^^ 



Z=0 
-1 



< cn 

It follows from (84) that the number of quantum queries satisfies 



n, 



(m-l \ 
n + J2 ( Z + l)2 íil '(n2-( dl+T )' + 1) 

< c(n + m2 dim ) < cn, 

and a simple change of variables in the index yields the desired result (46) 
for this case. 

Next assume min(s, d + a, d) = d\, put 



m 



logn 



(85) 



and 

N {0) = n, N 
Then we have, by (83), 



(1 = 0,..., m-l). 



e(T k ,A,£) 

< c(m ao 2- dim + m2- dm + n" 1 ) + 
m— 1 



1=0 



m—í 



< cmn _1 ^](l+l) 



1=0 



< cn 1 m 3 < cn 1 (logn)" 
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and, by (84), 



m— 1 



n, 



t(A) < c n+ ^(/ + l)2 cílí (nm- 1 2- dl ' + l) 



1=0 



< c(mn + m2 im ) < cnlogn. 

Again we get (46) by a change of variables in the index. 

Finally, we suppose min(s, d + a) < d±. Let r > be such that 



and put 



min(s, d + a) + t < di, 
logn 



m 



d 1 



7V(°)=n, Ni = \n2- dll - T(m - l ï] (l = 0, . . . , m - 1). 
We derive from (83) 
e(T k ,A,£) 

< c(m a "2~ ™in(s,d+<r)m + m2 --im + n ~l) + 

m— 1 

1=0 

< cn~ min (M+<x)M (bg n ) a ° + 

m—í 

c n -l2dil+T(m-l) ^ _|_ ^-ja 2-min(s,cí+o-)/ 

< cn^'^^'^^^logn) 00 

m— 1 

+cn _1 m ao 2 Tm 2( íil-T ~ min ( s ' c!+ÍT )) í 



< c^-minís^+^M ( logn )ao + cn -l m Q ° 2 (<*l -min( S ,d+ CT ))r 

< cn-^^'^/^Oogn)" . 
The number of queries is 



n, 



(m—l \ 
n + ( l + l)2 dll (n2- dll - T ( m ~ l ï + 1) 

< c^n^(/ + l)2- T ( m - í )+m2 c!im ^ < cnlogn, 
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and (46) follows. This concludes the proof of Proposition 7. □ 

Next we consider the case Qi = [0, l] dl and Q 2 = [0, l] d , where < di < 
d. Again we study the approximation of T k f with k G C s ' a (Q\, Q 2 ) a fixed 
kernel, but now / G C r (Q 2 ), and thus, the operator is considered as 
acting from C r {Q2) to L^Qi). Here rGN and C r {Q2) denotes the space 
of continuous complex- valued functions on Q which are r-times continuously 
differentiable in the interior Q®, and whose partial derivatives up to order r 
have continuous extensions to Q 2 . The norm on C r {Q2) is defined as 

ll/llc-(Q 2 ) = max sup \D a f(x)\. 

\ a \< r xeQ 2 

We let F = B C r (Q2) , G = Loc(Qi), S = T k , K = K, and 

A = {5 X : x G Q 2 }, (87) 

where <5 X (/) = f{x). 

To cover the logarithmic factors, we introduce for a G R with — d < a < 

+00 

if min(d + a,d) > di 
n(a) = < 4 if min(d + a, d) = d l (88) 
k if d + (7 < di = d, 

and if d + <r < di < d, we fix any eq > and define 



f 4 if < r -, + 1 



k(ct) 



i + 6 + £ if =^ = 5 + 1 (89) 
4 + £0 if r±|±£ >S + i. 



Proposition 8. Assume < di < d, s G N, s > di ; a G R, —d < cr < +oo ; 
rGN. T/ten í/iere is a constant c > suc/í í/iaí /or edZ G -Bc s . CT (Qi,Q 2 ) aní ^ 
n G N mí/i n > 2, 

eUT k ,B C r iQ2) ) < cn-H^^+i)^,)^. (90) 

where k((j) is as defined in (88), (89). 

Proof. Let n G N, n > 2. Let v be the smallest natural number such that 

e -i//8 < 1 / 8 _ ( 91 ) 

First we assume that either 

d + a > di (92) 
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(which, because of d\ < d, is equivalent to min(d + a, d) > d±) or 

d + a < di = d. (93) 

Comparing (88) with (45), we conclude that in these cases 

n(a) = P(a). (94) 

We write = T^J, with J the identical embedding C r (Q2) — ► C (Q2) , and 
Tk the operator T k , considered as acting from C(Q2) to L 00 (Qi). With 
X = C r (Q2), Y = C{Q2), and A as given by (87), the assumptions of 
Proposition 4 are easily verified. Therefore 

4un(Tk,B C r {Q2) )<8e^\J,B C r iQ2) ) sup e*(T k ,S). (95) 

í , e^(-B C (Q 2 )) 

By Proposition 7 and (94), 



sup el{f k ,£) < cn mm l^f - 1 ) (log n)^) (96) 
5G^(B C (Q 2 )) 

(the constants in this proof depend only on d,di,r,s,a). It is well-known 
that 

eT(J,B C r {Q2) )<cn- L *. (97) 
Furthermore, if (92) or (93) holds, we have 

r . fd+a \ . (r + d + a r \ 

d + mm V^T' V = mm V^^' d + 1 ) ■ 

This together with (95), (96), and (97) implies the desired result. 

Now we assume d + a < d\ < d (hence d\ / 0) . We recali the following 
construction from [18], proof of Proposition 4: We decompose Q2 = U£o 
with 

logn 



m 



d 1 

Hi = [0, l] dl x ^2-' [0, l] d - dl \ 2- (l+1 ï [0, l) d ~ dl ) (l = 0, . . . , m - 1) 

and 

H m = [0,l] dl x (y m [0,l] d ~ dl ) . 

Let 

k(x,y) = k(x,y) (x <E Qi, y <E Hi). 



(98) 
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Clearly, k e C s ' a {Qi,Hi) and 

\\ki\\c<°(Q 1 ,H l ) < \\k\\c»,°(Q u Q 2 )- (99) 
Put a\ = d\ — d. Arguing as in [18], proof of relation (52), we conclude 

Wh\\c^{Q,,H l )<c2^-^ 1 (0<Km). (100) 
We have the following representation 

m 

T k = Y,T h JiRu (101) 

1=0 

where Ri : C r {Q2) — ► C r (Hi) is the restriction operator, J\ : C r (Hi) — > 
C{Hi) is the embedding, and is considered as an operator from C{H{) to 
Loo(Qi)- With real numbers 5i,#2 > 0, which will be defined later, we put 



Vl 



} (^—Si)(m-l)-S 2 Í 



(0 < l < m). (102) 



Observe that p m = 1. Furthermore, define 

n , = (0 < / < m). (103) 

Note that ni > 2 for < l < m. As shown in [18], proof of Proposition 4, 
there is a constant c\ € N such that 

4tMu B C r {Hl) ) < c2- rl p- r (0 < l < m). (104) 

We verify that for < / < m 

sup el{T kv £)<c2~^ l pl d {\ogn l f. (105) 

£e^(Bc(H;)) 

Indeed, in the case < l < m relation (100) and Proposition 7 yield 

sup e%(T kl ,S) < cnf^^-^'Oogn,) 4 
5g^(-Bc(h ; )) 

< c2- dl Vr <Í 2 (<il_ ^ CT) '(logn í ) 4 

= c2-( d+,T ) z ^- cí (logn í ) 4 

If l = m, (99) and Proposition 7 give 



_ d+o- 



sup e£ m (T fcm , S) < cn m dl (log n m ) X' < c2~^ m (log n m ) X , 
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and, since p m = 1 and < 1, (105) follows. For l = 0, . . . , m we set 

vi = [8(21n(m-Z + l)+ln8)l . (106) 
It follows from (106) that 

m 1 m ^ 

£^ /8 %£(™-'+ir 2 <:r (107) 

1=0 1=0 

Define 

m 

fï = 2(ci + 2)i/ vin\ 
1=0 

with ci from (104) and v from (91). Then (101), (107) and Proposition 2 
imply 

m 

el( T k,B C r(Q 2 )) < 2 J2el {ci+2)uni (T kl JiRi,B C r {Q2) ). 

1=0 

The mapping Ri is of the form (1) with A = {5 X : x £ Q2}, A = {J^ : x £ 
íf/}, k = 1, and satisfies Ri{B C t^q 2 ^) Ç B C ri Hl \, hence, by Proposition 1, 

e 2( C1 +2)^n i ( T fci J '^' 5 C"·(Q 2 )) < e ? Cl +2) 1 ,n i ( T ^ J '' 5 C"·(// ! ))· 

Furthermore, by Proposition 4, 

Í +2)W1 C^J*,i?c W ) < Se^CJj.Bcr^,)) snp (T fcj ,£). 

Using this and (102-105), we get 

e n( T fc.- B C'·(Q2)) 
m 

< 16 E^Ji.BcW «up ^(T fe( ,£) 

z=o ££&{B c{Hl) ) 

m 

< cJ2^ (r+d+a)l P7 (r+d) (^gni) 4 

1=0 

'd. 



< cm 4 2~ (r+ (r+ )} (r+ 1 ) { ' (108) 

Relations (91), (98), (102), (103), and (106), give 

m m 

n = 2(ci + 2)vJ2"m < c^vfl 011 ^{m-l^dim-l^dl + l) , 
1=0 1=0 
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therefore, if 6\ > 0, 

- < c2 d im < m; ( 109 ) 

if ó 2 > 0, 

n < c2 dim \og{m + 1) < craloglog(n + 1), (110) 
and if <5i = 62 = 0, 

re < c2 dim mlog(m + 1) < cre lognloglog(re + 1). (HI) 

The proof will be accomplished by considering three cases. The first case is 
r + d + a < (r + d)^, that is, ,£± ^- < § + 1. Here we put S 2 = and take 
any <5i > satisfying 

r + d + a <(r + d)(^j-S^j . 

From (108) and (98), 

1 \ r-\-d-\-cr 

eUT k ,B cr{Q2) ) < cm * 2 -( r + d +^ m < cn-^(logn) 4 . 
Relation (109) and a suitable scaling lead to 

el{T k ,B C r [Q2) ) < crT^T^logn) 4 . 

The next case is r + d + a = (r + d)^j. Here we put S± = 62 = 0, and obtain 
from (108), 

eUT k ,B C r {Q2) ) < C m 5 2-^+ rf )^ m <cn-(5 +1 )(logn) 5 . 
Together with (111) this implies 

e^(T fc ,i? cr(Q2) ) < c (logn) 5 

vv ; \logreloglog(n + 1)/ 

< cn"(5 +1 )(logn)5 +6+£0 . 
Finally, if r + (i + <r > (r + d) ^ , we choose é>i = and 62 > so that 

r + d + <r - 82 (r + d) > (r + d) ^. 

From (108), 

el(T k ,B C r {Q2) ) < cmH-^ d ^ m < c^(5 +1 ) (log n) 4 , 
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4 



□ 



which together with (110) shows that 

*i(n,B C , m) < c( loglog ;,. +1) )~ (J+1 V g "> 

< cn-(5+ 1 )(logn) 4+£ °. 

6 Elliptic PDE 

Let d,m G N, d > 2, let Q C R d be a C°° domain (see, e.g., [18] for the 
definition), and let Jíf be an elliptic differential operator of order 2m on Q, 
that is 

j^u= Yj aa(x)D a u(x), (112) 

\a\<2m 

with boundary operators 

38ju= ^2 b ja{x)D a u{x), (113) 

|a|<mj 

where j = 1, . . . ,m, rrij < 2m — 1 and a a G C°°(Q) and ò JQ , G C°°(dQ) 
are complex-valued infinitely differentiable functions. We study the homo 
geneous boundary value problem 

Sfu(x) = f(x) (x G Q°) (114) 

SSju{x) = (x G <9Q). (115) 

Let 

a(x,£) := a a (x)£ Q (x G Q, Ç G M d ) 

|a|=2m 

b^x.O := b ja (x)C (xedQ, ÇeR d , j = l,...,m). 

\a\=rrij 

We assume the ellipticity condition: 

a(x,£)^0 (x G Q, (£1^ {0}) 

and for all linearly independent £, rj G M d the polynomial a(x, £ + rry) has 
exactly m roots r t + (i = 1, . . . , m) with positive imaginary part. Put 

m 

a + (x,Ç, V ,T)=Y[(T-T+). 

i=l 
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We also assume the complementarity condition: For all x G dQ and all 
Cxi v x G ^ d \ {0}, where Ç x is tangent to at x and is orthogonal to the 
tangent hyperplane at x, the set of polynomials bj(x, £, x +tv x ) (J = 1, ■ ■ ■ , m) 
is linearly independent modulo a + (x, Ç x , v x , r). Finally we suppose that 
there is a kq with < k < 1 such that for all / in the Hòlder space C K ° (Q) 
the classical solution u exists and is unique (see [20] and also, e.g., [1], for 
the assumptions made here). 

Let M be a smooth submanifold of Q of dimension di, where < d\ < d 
(see, again, [18] for a definition). If d\ = 0, we assume M = {x}, where 
x is any inner point of Q. Let r G N, F = B C r^, G = L 0Q (M), and let 
5 : F — > G be given as 

Sf = u\ M , (116) 

where n is the solution of (114), (115). So we want to find an approxima- 
tion of the solution of (114), (115) on a tii-dimensional submanifold M of 
the domain Q, for right-hand sides belonging to B C r^, and the error is 
measured in the L oc (M) norm. We put K = C and 

A = {8% : x€Q,\a\<r}, (117) 

where 5 x {f) = D a f(x), that is, we allow information consisting of vàlues of 
/ and its derivatives up to order r. 

Theorem 1. There are constants ci,C2 > such that for all n G N with 
n>2, 



cin v d i d / < e^(5, F) 

< ^-^(^'^(logn)^- d \ (118) 
wií/i k as defined in (88) and (89). 

Proof. By a result of Krasovskij [20], Theorem 3.3 and Corollary, there is 
a kernel k G C°°' 2m - d {Q,Q) such that for all / G C K °{Q) the solution u of 
(114), (115) satisfies 

u(x)= [ k(x,y)f(y)dy (x G Q). (119) 

Consequently 

(Sf)(x) = (T k f)(x) (x G M), 



37 



that means, S = Tk, with Tk considered as an operator from C r {Q) to 
L oa (M). First we prové the upper bound. We show that it holds even for the 
smaller sets of information functionals A = {ó x : x G Q}. Let Q\ = [0, l] dl , 
considered as a subset of M. d by identifying IR^ 1 with M. dl x {0( d ~ dl ^}, and 
let Q2 = [0, l] d . The following representation of 7\ was shown in [18], proof 
of the upper bound in Theorem 1: There is a p G N (depending only on M 
and Q) such that 

p 

Tk = Y, (Xi T ki Y i + XiTfH J) • (120) 
i=i 

Hcrc 

Xt : Loo(Qi) - ^oo(M) 
: C r (Q) - C r (Q 2 ) 

(i = are bounded linear operators and J : C r (Q) — ► C(Q) is 

the embedding. Moreover, 1^ is of the form (1) with A = {5 X : x G Q}, 
A = {5 X : x G Q2}, and re = 1. The kernels satisfy 

h G C 00 ' 2m - d (Q 1 ,Q 2 ), (121) 
hi G f| C^Q^Q), (122) 

cr>0 

the integral operator is considered as acting from C r (Q2) to L OCJ (Qi), 
and is considered as a mapping from C(Q) to L oc (Qi). (Using the 
terminology of [18]: up to shifting and scaling of cubes, Xj stands for the 
product EíXí, hi is the k[ from relation (64) of [18], and h{ is k" from relation 
(66) of that paper, extended by zero to all of Q.) 

Let z/o, z/i be the smallest natural numbers satisfying 

pe -V8 < 1 /g ) e -^i/8 < 1 / gj 

respectively. Let Ci = z/o(3z/i + 2)p. Then, by Proposition 2, 

e cm( T k, B Cr(Q)) 
p 

< 2 2 (eUX t T k Y t , B C r (Q) ) + e^J^T^J ', £ C . (Q) )) . 
í=i 

By Proposition 1, 

e 2n( X i T ki Y i, B Cr(Q)) < \\Xi\\e%(T k ., ll^ll-^C r (Q 2 )) 

= WXilWlYilKiT^Bcr^). 
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Furthermore, by Proposition 4, 

e L· 1 n(XiTh i J,B C r(Q)) 

< 8e^(J,B Cr(Q) ) sup ^(XiTh.,8). 

£6^(Bc(Q)) 

Moreover, Lemma 1 of [14] gives 

sup e%(XiT hi ,£) < \\Xi\\ sup e%(T hi ,£). 

£e c <f(B C (Q)) £^( B c(Q)) 

Thus we obtain 

e c in ( T k,B C r {Q) ) 

< cY,{el(T ki ,B Cr{Q2) ) 

i=i 

+e* rt (J,B Cr(Q) ) sup e<*(T hi ,£)). (123) 

We conclude from (121) and Proposition 8 that 

eï(T ki ,B C r {Q2) ) < c^ min (T'5 +1 )(logn) K ( 2m - d ), (124) 
where /t(2m — <i) is defined in (88). Moreover, 

e d n et (J,B C r {Q) )<cn-i, (125) 
see, e.g., [29]. Furthermore, from (122) and Proposition 7, 

sup e*(T hi ,£) Kcn-^logn)^ (126) 



where 



, if di < d, . , 

ai = ^ 4 if dl = d. (127) 



Relations (123)-(127) finally give 

r+2m r _ 



e ^ in (r fc ,5 c , (Q) ) < cn" mÍ H^r-^ +1 J(logn) K ( 2m - d ). (128) 

Indeed, in the case d\ < d this is clear. For di = d we argue as follows: If 
d + (2m — d) = 2m > d = d\, then by (88), n(2m — d) = 4, which gives 
(128). If 2m < d = d u then = r±*m < § + 1, hence 



n 



"(d +1 )(logn) 4 < cn 
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which leads to (128), again. Now the desired upper bound in (118) follows 
from (128) by rescaling. 

Next we prové the lower bound. As above, let Q\ = [0, ï\ dl , Q 2 = [0, ï] d . 
Let Cq(Q 2 ) denote the subspace of C r {Q 2 )) consisting of those functions 
whose partial derivatives up to the order r vanish on the boundary of Q 2 . 
It was shown in [18], section 5, that there are bounded linear operators 
*o : C T {Q 2 ) -► C r {Q) and Y : L^M) -► C such that 

Y SX = S u (129) 

with Si : Cq(Q 2 ) — > C the integration operator 

Sif= í f(y)dy (feC r (Q 2 )), 

and Xo is of the form (1) with A = {6% : x G Q 2 , \ct\ < r}, A = {5% : x G 
Q, | ck | < r}, and k depending only on d and r. Consequently, by Proposition 
1 and (129), 

4^1,^2)) < ll^o||e^S,||X ||^ {Q) ) 
= \\X \\\\Y ü \\el{S,B C r {Q) ). 

From [23] it is known that 
Thus we conclude 

el(S,B Cr(Q) ) > cnS- 1 . 

This proves the lower bound of (118) for the case ^ + 1 < r+ ( ^ m (including 
the case d\ = 0). 

Now we assume d\ > 1. We use another reduction from [18], section 5, 
which will give the remaining part of the lower bound: There are bounded 
linear operators X : Cg +2ín ((5i) — ► C r (Q) (representing the composition 
££XE from [18]) and Y : L^M) -► Loo(Qi) such that 

YSX = J, (130) 

where J : C +2m (Qi) — > L OQ (Qi) is the identical embedding, and X is of 
the form (1) with A = {5" : x G Q®, a G Ng\ \a\ < r + 2m}, À = {££ : 
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x G Q, a G Nq, |q| < r}, and k depending only on d\, m and r. From 
Proposition 1 we obtain 

e* Kn {YSX,B c[ + 2m(Qi) ) < \\Y\\eUS,\\X\\Bcr( Q) ) < ce*{S, B Cr{Q) ) . 

Together with (130) this yields, 

e 1 K n( J > B c{+ 2m (Qi)) ^ ce n{S,B C r {Q) ). 



By [15], 
Consequently, 



<= q ( T fí j_9 ï > r í7 -( r + 2m )/' i 
e 2 K nW'- tí C£ +2m (Qiy - Cn 



eUS,B Cr(Q) )>cn-^ 2m y d \ 
concluding the proof of the lower bounds. □ 

7 Comments 

In this section we recali previous results on the complexity of elliptic equa- 
tions in the classical deterministic and randomized setting and compare them 
with the results of the present paper. We discuss the speedups between the 
different settings. 

Below S and F refer to the elliptic problem studied in section 6, see 
(116). Let e^ et (5, F) and e^ n (S, F) be the n-th minimal deterministic and 
randomized errors, respectively, as introduced, e.g., in [18], section 3. 

To suppress logarithmic factors, we use the following notation: For func- 
tions a, b : No — > [0, oo) we write a(n) >q g b(n) if there are constants 
ci, C2 > 0, riQ e No, «i, «2 G M such that 

ci(log(n + 2)) Ql 6(n) < a{n) < c 2 (log(n + 2)) a ' 2 b(n) 

for all n G No with n > uq. Furthermore, we write a{n) x b(n) if the above 
holds with a.\ = «2 = 0. 

In the classical deterministic setting we have 

e f(S,F)xn-Í 

This result is essentially contained in [30, 31, 5], see also [18], where a proof 
is given for the specific function spaces considered here. Observe that in the 
deterministic setting the rate does not depend on d\, thus local and global 
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problem are (up to constants) equally difficult. As established in [18], in the 
classical randomized setting we have 



e™(S,F) x log n 
By Theorem 1, in the quantum setting, 

e«(S,F) x log n 

Thus, as in the classical randomized setting, the rate in the quantum setting 
depends on d\. Note that the rate is n _ d _1 (the same as that of quantum 
integration of functions from C r (Q), see [23]), for all d\ < min(eZ, 2m) and 
r G N. Indeed, if d < 2m, we infer 

r + 2m r + 2m r 

— ■ > — ■ > - + 1, 

di d d 

while, if d > 2m and d± < 2m, we have 

r + 2m ^ r + 2m r 
d\ ~ 2m d 

Let us compare the quantum setting with the classical deterministic setting. 
We have a speedup for all < d\ < d: For example, for all d\ < min(d, 2m), 
the speedup is n _1 . Furthermore, if d\ = d and d > 2m, the speedup is still 

2m 

n <i . 

Comparing the quantum with the classical randomized setting, we see 
that for di < min(d, 2m) there is a speedup of n~ï, while for d\ = d and 
2m < d < 4m the speedup is nï 5", and there is no speedup at all for 
di = d and d > 4m. 
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